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Abstract 


The extreme rate of evaporation of the Dead Sea (DS) has serious implications for 
the surrounding area, including atmospheric conditions. This study analyzes the aerosol 
properties over the western and eastern parts of the DS during the year 2013, using 
MAIAC (Multi-Angle Implementation of Atmospheric Correction) for MODIS, which 
retrieves aerosol optical depth (AOD) data at a resolution of 1km. The main goal of the 
study is to evaluate MAIAC over the study area and determine, for the first time, the 
prevailing aerosol spatial patterns. First, the MAIAC-derived AOD data was compared 
with data from three nearby AERONET sites (Nes Ziona - an urban site, and Sede Boker 
and Masada - two arid sites), and with the conventional Dark Target (DT) and Deep Blue 
(DB) retrievals for the same days and locations, on a monthly basis throughout 2013. For 
the urban site, the correlation coefficient (r) for DT/DB products showed better 
performance than MAIAC (r=0.80, 0.75, and 0.64 respectively) year-round. However, in 
the arid zones, MAIAC showed better correspondence to AERONET sites than the 
conventional retrievals (r=0.58—0.60 and 0.48-0.50 respectively). We investigated the 
difference in AOD levels, and its variability, between the Dead Sea coasts on a seasonal 
basis and calculated monthly/seasonal AOD averages for presenting AOD patterns over 
arid zones. Thus, we demonstrated that aerosol concentrations show a strong preference 
for the western coast, particularly during the summer season. This preference, is most 
likely a result of local anthropogenic emissions combined with the typical seasonal 


synoptic conditions, the Mediterranean Sea breeze, and the region complex topography. 
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Our results also indicate that a large industrial zone showed higher AOD levels compared 
to an adjacent reference-site, i.e., 13% during the winter season. 
Graphical Abstract 


Different AOD patterns emerged over the Dead Sea coasts. 


Keywords: Aerosols, Aerosol Optical Depth (AOD), MAIAC, MODIS, AERONET, the 
Dead Sea. 


1. Introduction 


The Dead Sea region is the lowest place on earth, located at a height of 
approximately 430 meters below sea level (BSL), and has unique climatic and landscape 
characteristics. Over the last 2,000 years, the region has changed dramatically (Bookman 
et al., 2004; Kottmeier et al., 2012); suffering from an extreme rate of evaporation, the 
Dead Sea is drying out, and its surface level has dropped some 40m in the last 40 years 
alone (Alpert et al., 1997; Shafir and Alpert, 2011). Furthermore, the Dead Sea area is 
characterized by a near-constant layer of haze owing to the extensive evaporation, which 
in turns leads to a more obstructed pathway for incoming solar radiation (Even-Paz and 
Shani, 1989). These conditions mean that the area is an important one for studying 
extreme aerosol pollution, which can reach a concentration of up to several thousand 
micrograms per cubic meter and have a negative impact on human health (Kishcha et al., 
2016). On the extreme dust event of September 2015, concentrations reached a maximum 
of 1700-10,000 ug m-3 (based on three monitoring stations on the western side) (Alpert 
et al., 2016, Gasch et al., 2017; Mamouri et al., 2016), whereby the average annual value 


is around 50 ug m-3 (Kishcha et al., 2016). 
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The following factors influence aerosol dynamics in this region: 1) the very hot 
and dry climatic conditions (Shafir and Alpert, 2011); 2) large quantities of aerosols from 
different sources; 3) continual decline in the water level of the Dead Sea, which weakens 
the local Dead Sea breeze while intensifying the Mediterranean Sea breeze penetrating 
the valley (Shafir and Alpert, 2011); and 4) natural evapotranspiration of salts, bromine, 
and other minerals. Therefore, the study of aerosols in the Dead Sea can greatly enhance 


our understanding of extreme aerosol pollution and our ability to predict it. 


Previous studies on the region’s haze and dust aerosols have focused on a wide 
range of topics. Moore et al. (2013) studied spatial and temporal patterns of atmospheric 
mercury depletion events. Singer et al. both examined dust deposition using collectors 
installed at ground level (Singer et al., 2003), and the properties of suspended dust by 
means of direct sampling (Singer et al., 2004). Suspended atmospheric particles were also 
characterized by Kalderon-Asael et al. (2009). The vertical distribution of haze particles 
with regard to their physical and chemical properties was described by Levin et al. 
(2005). With other studies examining the relations between mineral dust aerosols and air 
quality (Ganor et al., 2009), the area’s atmospheric optical depth and UVB irradiance 
components (Kudish and Evseev, 2012), and the dust transport in the region by modeling 
(Vogel et al., 2006). However, to the best of our knowledge, no research has ever been 
conducted over the Dead Sea coasts to identify variability in aerosol levels. Furthermore, 


efforts to study this variability are hindered by the lack of ground monitoring data. 


Satellite imagery is an important tool for evaluating air quality and pollution, 
owing to the extensive spatio-temporal monitoring of the earth surface and atmosphere 


(Kaufman et al., 2002; Kokhanovsky et al., 2007; Kaskaoutis et al., 2012). The main 


parameter derived from satellite observations for assessing air quality and atmospheric 
opacity is aerosol optical depth (AOD), a measure of the extinction of electromagnetic 
radiation at a given wavelength due to the presence of aerosols in an atmospheric column 


(Kaufman et al., 2002). 


Since the Dead Sea is located in an arid zone, the bright arid areas surrounding it 
constitute a challenge for using remote sensing data due to the difficulty in discerning 
aerosols above the bright background terrains (Remer et al., 2005). With the release of 
the “Deep Blue” (DB) product, which is based on data from the 412nm, 470nm, and 
670nm spectral channels, it has become possible to study the spatial and temporal 
variability of transported dust/pollution over deserts, as bright surfaces show low 
reflectivity in the blue spectral region (Hsu et al., 2013). These data are available from 
the Sea-viewing Wide Field-of-view Sensor (SeaWiFS) (in operation from September 
1997 to December 2010) and the Moderate Resolution Imaging Spectroradiometer 
(MODIS) (2000-current), which employ the DB algorithm (Hsu et al., 2004; Hsu et al., 
2006). Furthermore, the recent MODIS Collection 6 (C6) aerosol retrievals include 
enhanced 10km DT and DB AOD and a “merged” DB-DT parameter at 10km and 3km 


resolutions (Levy et al., 2011, 2013; Hsu et al., 2013). 


A Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm 
(Lyapustin et al., 2011b) for MODIS provides a 1km resolution aerosol retrieval. MAIAC 
has been shown to be successful in obtaining detailed AOD data over bright urban areas - 
in New England (Lee et al. 2011; Chudnovsky et al, 2013a; Chudnovsky et al. 2014), the 
southern United States (Lee et al., 2016), across the Mexico City metropolitan area (Just 


et al., 2015), and several regions in Israel (Kloog et al., 2015). MATAC also has the novel 
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ability to detect biomass burning (smoke) aerosol (Lyapustin et al., 2011b). The 
algorithm was tested in retrieving AOD in the European Alpine region, and showed good 
correlations (for MAIAC filtered data) with AERONET (AERosol Robotic NETwork) 
measurements (R ranging from 0.78 to 0.90) (Emili et al., 2011). Another advantage of 
MAIAC is its improved detection of cloudy and clear-sky areas, enabling it to produce 
AOD data over partially cloudy days, when the conventional 10km resolution MODIS 


AOD returns no retrievals (Chudnovsky et al., 2013a, 2013b). 


1.1 Problem definition: The need for high-resolution aerosol monitoring 


With such extreme and rapid changes affecting the Dead Sea area, the 
environmental monitoring of aerosols is very important. As clearly seen in Figure 1, the 
conventional MODIS aerosol product (a 10km resolution) is insufficient for identifying 
aerosol patterns in this area. A solution may be offered by MAIAC (Panel B), which is 
able to retrieve data in this complex and ever-changing terrain, making it possible to 


identify variability in AOD patterns. 


[INSERT FIGURE 1 HERE] 


Given the challenges described above, the main goal of our study was to assess 
the accuracy of MAIAC over the Dead Sea coasts and, for the first time, determine the 
concentration patterns of aerosols in the region. To that end, the association between 
monthly data (January-December 2013) from MAIAC-derived AOD and from the three 
AERONET sites was investigated, and correlations to the conventional DT and DB 


retrievals for the same days and locations were calculated. In addition, we generated 
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seasonal AOD maps and conducted a year-long study of the differences between AOD 
levels on the eastern and western coasts of the Dead Sea. This enabled us to identify 
distinct aerosol patterns over the Dead Sea coastal region. Finally, we compared 
differences in AOD levels between one of largest anthropogenic sources of pollution 


(Rotem factory) and a nearby reference site (each area consisting of 14 pixels). 


2. Study Area 


The area chosen for this study comprises the land area immediately surrounding 
the Dead Sea. (The study of aerosols above the sea itself was excluded from the analysis 
due to an insufficient number of retrievals over the water). Hence, an area of roughly 


120km x 150km was analyzed (Figure 2). 


[INSERT FIGURE 2 HERE] 


In addition to the strip of hotels along the southern part of the sea, there are also 
several village/rural settlements and medium-sized cities in the region. The Dead Sea 
Works (potash plant) is located on the southwestern coast of the Dead Sea (at 390m 
BSL), and the Rotem factory (400m ASL), belonging to a large mining company that 
extracts phosphates from the Dead Sea and manufactures agricultural fertilizers and other 


chemical products, is found further to the west. 


3. Material and Methods 


3.1 Satellite data 
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Our study relies on MAIAC AOD retrieval (at 470nm) from the Aqua satellite 
MODIS (13:30 local time). MAIAC is an algorithm developed for MODIS that produces 
AOD data at a 1km horizontal resolution (regardless of scanning angle) by using a time 
series of measurements from MODIS while simultaneously processing a group of pixels, 
with regional background models based on AERONET climatology (Lyapustin et al., 
2011a). The aerosol retrieval procedure is based on the spectral similarity of the surface 
Bidirectional Reflectance Distribution Function (BRDF) shape between the visible 
(0.47um) and shortwave IR (2.13um) wavelengths, with look-up tables that contain pre- 
computed radiative transfer functions (Lyapustin et al., 201la, 2011b, and 2012). 
Detailed information regarding MAIAC’s radiative transfer basis and the look-up tables, 
aerosol algorithm, and atmospheric correction is described extensively by Lyapustin et al. 
(201 1a, 2011b, and 2012) in a number of publications. 

Cross-validation of MAIAC data with AERONET has showed MAIAC’s superior 
accuracy in retrieving AOD data over bright (urban) surfaces in comparison with the 
standard MODIS AOD product (dark target algorithm), stemming from its more detailed 
surface characterization method (using surface BRDF rather than Lambertian reflectance 
model) and cloud detection (Lyapustin et al., 2011b). On the other hand, individual 1km 
retrievals may sometimes have added noise, due to residual uncertainties in the surface 
reflectance or from cloud/snow contamination (Emili et al., 2011, Lyapustin et al., 201 1a, 


2011b, and 2012). 


In addition to MAIAC data, we used the Collection 6 combined aerosol product 
that includes Dark Target (DT) and Deep Blue (DB) AOD data for the period of 2013 at a 


spatial resolution of a 10kmx10km (at nadir). This product, over land and water, aims to 


fill the gaps in the three individual algorithms, creating a more complete picture of 
aerosols at 10x10km resolution and at 550nm (Sayer et al., 2014). Aqua MODIS C6 
includes an operational combined AOD product calculated from DB and DT AOD in 
three Normalized Difference Vegetation Index (NDVJ) categories (Hsu et al., 2006). This 
combined AOD is equal to DT AOD if NDVI > 0.3 and is equal to DB AOD if NDVI < 
0.2. Detailed descriptions of the MODIS operational combined AOD algorithm and the 


QA flags can be found in Levy et al. (2013). 


3.2 AERONET observations 


Our analyses were geographically extended to include all available AERONET 
stations with high temporal resolution located in Israel (Figure 1): the Nes Ziona station, 
located in an urban region (part of the Tel Aviv metropolis) at a height of about 40m 
(ASL); the Sde Boker station, located in a kibbutz in the northern part of the Negev 
Region (an arid zone) at a height of 480m ASL (this station is the closest one to the Dead 
Sea); and the Eilat station, located at the southernmost point of Israel, near the sea shore 


of Aqaba bay (also an arid zone) at a height of (15m ASL). 


The AOD values from the AERONET stations were chosen for the times 
correlating to the satellite passing over the region. Since not all AERONET AOD 
measurements were performed during the satellite pass time, eventually all analysis was 
conducted using the average AOD value of a one-hour time interval, from half an hour 
before the satellite pass to half an hour afterwards (13:30 ascending for Aqua). For each 
station, the AOD distribution was matched to the MAIAC AOD distribution, comparing 
the AOD value at the AERONET point location (a given latitude and longitude) to the 


AOD value of the MATAC lkm x 1lkm pixel where the AERONET station was located. 
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With regard to wavelength consistency, MAIAC provides AOT at 470nm, where the 
retrievals are made, and at 550nm for the users' convenience. In this work, AERONET 
AOT at 550nm was obtained by linear interpolation which has a sufficient accuracy for 
this study. A higher accuracy could be obtained using second order interpolation (e.g., 


Eck et al., 1999). 


3.3 Data processing and analysis: the methodological approach 


Data processing was conducted in four major stages (Figure 3). First, we have 
constructed a data base for analysis for the year 2013. All the retrievals from MAIAC 
were subjected to initial screening for days suitable for our study. AOD data for the three 
AERONET station located in Israel and for MODIS standard product was downloaded as 
well. All three data sources were combined to one data set by date. Second, in order to 
validate MAIACs quality of retrieval in the region, we have done a direct comparison of 
AERONET data to satellite-derived DT/DB and MAJAC data. Third, in order to identify 
seasonal differences in AOD concentration and patterns we have generated seasonal 
AOD maps. Next, we have compared aerosol pollution levels between the eastern and 
western areas of the Dead Sea. Lastly, we compared AOD levels between the large 


industrial zone of Rotem factory to a nearby reference site. 


[INSERT FIGURE 3 HERE] 


3.3.1 Data screening 
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All AOD retrieved images collected during 2013 were examined and compared to 
RGB images in order to identify cloudy days and days with unexplained missing data. 
This visual comparison was done since MAIAC have never been tested in the 
topographically complex region of the Dead Sea. To that end, a program in MATLAB 
language that opens and plots AOD and RGB images for each day simultaneously was 
written. This manual daily analyses for relatively new algorithms was also done by Ben- 
Ami et al. (2009) when CALIOP data was firstly released. Figure 1s shows examples of 
days that were excluded from our analysis, and days that were eventually included in the 
analyses despite representing some inaccuracies. The most common reasons for 
disqualifying data (17% of the sample in total was disqualified) were cloud cover (as in 
Panel B, for August 30, 2012) or missing AOD retrievals (as in Panel A, for August 3, 
2012). In Panel A, we see clearly an area bordering the Dead Sea to the east (marked by a 
white circle) for which no data was retrieved and where no significant clouds were 


identified in the RGB image. 


Data that exhibited residual calibration errors in the original MODIS L1B data, 
related to polarization effects (Franz et al., 2008), were included in our analyses. This 
phenomenon can be identified by parallel semi-horizontal stripes across the image. 
Sometimes these lines cannot be fully removed, even after a polarization correction is 


performed, as is demonstrated in Panel C of Figure 1s (August 8, 2013). 


3.3.2 AERONET vs satellite-derived AOD 


Next, we investigated the associations between satellite-AOD and ground-based 


AERONET measurements (see Figure 2 for AERONET station locations) for the year 


11 


2013. We made a comparison between Collection 6 DT, DBMYD04, and MAIAC 
retrievals for all available days on a seasonal basis. Following, we restricted our analyses 
to days when all retrievals were available for a given AERONET site (that is, same days 
and locations). Specifically, for the urban area we compared MAIAC, DB, and DT, 


whereas for the arid zone we compared MAIAC and DB. 


In addition, we examined whether the relation between ground-based AERONET 
measurements and satellite AOD retrievals becomes stronger as the spatial resolution of 
the AOD becomes coarser. To that end, we used different spatial resolutions to determine 
what “box size” (representing a given number of pixels) would improve the correlation 
between MAIAC and AERONET measurements. The analysis was repeated for 
progressively degraded resolutions, from a 3km x 3km box (9km7), to a 5km x 5km box 
(25km’), and a 10km x 10km box (100km’), obtained from the original 1km AOD data 
by simple averaging. Here we assumed that the coarser resolution can be adequately 


represented by a simple aggregation of the MAIAC fine resolution retrievals. 


3.3.3 Seasonal spatial pattern and differences in aerosol levels between Dead Sea 


coasts 


Following data screening, seasonal AOD averages were produced for all days that passed 
through our data quality test. Next, to capture the difference between two locations i and 
k, referred to as the Normalized AOD Difference (NAODD), we divided the difference in 
AOD values between both locations for the same days by the sum of their corresponding 


AODs, as follows: 
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NAODD= (AOD , ;- AOD ,,) / (AOD ,;+ AOD ,,) (1) 


where j is the date of observation (that is, each point represents the calculated normalized 


difference on a given day). 


Negative values of NAODD would correspond to days when location k is more 
polluted than location i, while values close to zero would correspond to days when 
pollution levels at locations i and k are similar, and high positive values would indicate 
days when location i is more polluted. Here, we compare between the eastern and western 
coasts of the Dead Sea. In addition, we compared between the Rotem factory site to a 
nearby reference location located at a height of 400 meter above sea level. To that end, 
monthly average AOD was calculated for both locations. Here, only months with highest 


difference in average AOD values are reported. 


4. Results 


4.1 Direct comparison between DT/DB, MAIAC and AERONET retrievals 


This section studies the subset of MYD04/MAIAC data for days when retrievals 
were available for a given AERONET site. Table 1 shows the comparison of AOD from 
AERONET with AOD from (MODIS) DT/DB and MAIAC in Israel (three locations) for 
2013. As can be seen, the correlation varies by site. For the urban site, MAIAC shows 
slightly lower yearly correlations than DB/DT retrievals. Note that if dust storm days or 


“extreme points” were included, the correlations would have been similar and higher. 


[INSERT TABLE 1 HERE] 
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Table 2 shows a direct comparison between AOD from AERONET with DT/DB 
and MAIAC on a seasonal basis, for the same days in 2013 and the same three locations 
in Israel. While DT performs better in summer for the urban site, MAIAC is more 
accurate for the arid zone, producing better correlations there for all seasons except 


winter in Eilat. 
[INSERT TABLE 2 HERE] 
4.2 AOD pattern 


Figure 4 shows monthly MAIAC and AERONET time-series measurements for 
the Sde Boker site during four representative months: April, July, August, and November 


2013. 
[INSERT FIGURE 4 HERE] 


As can be seen, there are days when MAIAC overestimates or underestimates 
AOD compared to ground-based observations. However, the general monthly trend and 
variability of AOD from satellite-derived MAIAC are quite similar to those from 
AERONET. With this result in mind, we calculated seasonal means of AOD retrievals; 


these are presented in Figure 5. 
[INSERT FIGURE 5 HERE] 


Figure 5 shows the seasonal averages of MAIAC AOD for 2013. During spring, 
the dustiest season, there are higher AOD values over the entire region. The winter 
(December—February) and fall (September-November) months have a lower seasonal 
mean AOD, while during the summer we see a concentration of somewhat higher values 


around the western coast of the Dead Sea, where the northwestern wind penetrates the 
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Negev, bringing humidity and air masses from the densely-populated industrial coast of 
central Israel as well as from Europe (Derimian et al., 2006). Previous studies has found 
that during this season and at this location, fine particle loadings with anthropogenic 


origin are dominant (Derimian et al., 2006; Dayan et al., 2017). 


As shown in Figure 1, during selected days in August, MAIAC exhibits high 
capability in retrieving AOD data over the bright arid area of the Dead Sea. This is 
apparent from the higher spatial coverage when comparing the algorithm to the 
conventional MODIS product. To further investigate the variability in AOD patterns over 
the region and their seasonal dependence, we calculated the NAODD (Normalized AOD 
Difference index) for 2013, and the results are shown in Table 3. As our results indicate, 
the western and eastern parts of the Dead Sea area have slightly different trends in 
NAODD values, lightly skewed toward positive values (that is, the west coast is more 
polluted). The 25th percentile is -0.22, the 50th percentile is 0.065, and the 75th 


percentile is 0.16, with a mean NAODD of 0.074. 


[INSERT TABLE 3 HERE] 


On average, the western part of the Dead Sea has consistently higher aerosol 
concentration than the eastern side (Table 3). During the spring season (67 available days 
for analysis), more than half of the days in the sample (56.7%) exhibited higher AOD 
levels on the western coast. This percentage increased during the fall season (a sample of 
77 days), where 62.3% of the days had higher AOD levels on the western side than on the 
east. This trend is seen at best during the summer season (91 days), when on 92.3% of the 
days, the aerosol-pollution levels were higher on the western side. During the winter, 


both coasts show similar AOD levels, with pollution being higher on the eastern/western 
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side on 53.3/46.7% of the days, respectively. These results are not surprising and can be 
explained by the prevailing synoptic conditions in the region. In the summer, the 
Mediterranean Sea breeze is very strong on the western side and significantly weakens on 
the eastern side (Shafir and Alpert, 2011), which can explain the pronounced differences 
between the east and the west. In contrast, in the winter, the strong prevailing winds and 
the activity of the winter lows may create the observed similarity between the two coasts. 


While in the fall and spring intermediate synoptic conditions prevail. 


Simulating the wind regime using the WRF model (three-dimensional scale), we found 
further explanation to our results (personal communication confirmed by numerical 
simulations, with Kunin, P., Rostkier-Edelstein, D., and Alpert, P.). The simulated local 
dynamics show the entrance of the Mediterranean Sea breeze at the western side, with a 
strong downward vertical velocity followed by a strong upward vertical wind, pertaining 
to the area topography. However, on the eastern side the upward motion is fairly weak. 
This wind dynamics was also discussed in the work of Bitan (1982), where he analyzed 
10 years of data from 15 meteorological stations, and clearly identified the strong 
penetration of Mediterranean Sea breeze in the west during the summer (Figure 2, in 
Bitan, 1982). Therefore, the wind regime described above bring about the strong local 
dust and aerosols emissions, with seasonal dependency which creates the pronounced 
differences between winter and summer. Several examples of days with variability in 


AOD levels between the two coasts are shown in Figure 6. 


[INSERT FIGURE 6 HERE] 


Would a coarser resolution yield a better proxy for the AERONET observations 


measured at the ground? Here we assumed that perhaps the surface reflectance might 
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have an impact on the accuracy of MAIAC AOD retrievals. The answer to our question is 
presented in Table Is. As can be seen, the change in resolution does not affect the 
accuracy of urban Nes Ziona site, but does slightly improve the correlation for Sde Boker 
site. This result indicates that a coarser resolution contribute very little if at all to the 
improvement in the correlation between MAIAC and AERONET AOD measurements. 
Therefore we can not recommend using this tool as a standard procedure when 
anaylysing AERONET data, and its impact should be examined by site. 

Finally, we compared between Rotem factory to a reference area with similar 
topographic conditions. Our results indicate that this large industrial zone showed higher 
AOD levels compared to an adjacent reference-site, i.e., 13% during the winter season. 
Other periods exhibited much lower differences (April-May, 2-5% correspondingly) or 


even the opposite result (March). 


4. Discussion 


In this study, we examined the utility of MAIAC retrieval to assess the aerosols 
vertically-integrated concentration (AOD) over the arid region of the Dead Sea and the 
prevailing aerosol patterns. As seen in our presented results, the AOD patterns (Figures 1 
and 6) observed for the Dead Sea area not only indicate to the strength of MAIAC in 
retrieving data above bright reflecting regions, but also suggest that the region generates 
typical re-occurring patterns. In addition, there is a relatively good agreement between 
the MAIAC retrievals and AERONET measurements, as apparent from Figure 4 for the 


nearby Sde Boker AERONET site. 
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The aerosol patterns seen over the Dead Sea region mostly show higher 
concentrations over the western part of the area compared to the eastern part, in particular 
during the summer season. This result is not surprising and can largely be explained by 
four main factors: (1) the unique synoptic and meteorological conditions; (2) the 
Mediterranean Sea breeze effect (Bitan et al. 1982; Kishcha et al., 2016; Alpert et al., 
1997); (3) the breeze created by the Dead Sea itself; and (4) the areas unique topography 


and location (Bitan, 1974, and Bitan, 1977), as follows. 


During the summer season, the persistent Persian Trough system cause local 
stability beneath the seasonal strong inversion layer at about | km above mean-sea-level. 
This allows the Mediterranean sea breeze to influence the area greatly alongside its 
unique topography, by enhancing the impact of local emission sources (Kishcha et al., 
2016; Alpert et al., 1997; and Dayan et al., 2017) and transporting dust and 
anthropogenic aerosols from afar (Derimian et al., 2006; Derimian et al. 2017). 
Specifically, Derimian et al. (2017) have shown that chemical composition, 
microphysical and optical properties of atmospheric aerosol deep inland in the Negev 
Desert of Israel were found to be influenced by daily occurrences of sea breeze flow from 
the Mediterranean Sea. A similar effect was described by Tokar et al. (1993), where 
authors demonstrated how pollution in Hadera (a coastal city in Israel) progresses inland 
due to this wind regime. During winter, however, strong winds and winter lows activity, 
contribute to stronger mixing and therefore yield to similar AOD levels on both coasts. It 
is interesting to note that a similar process was described for the region of lake Kinneret 
(the Sea of Galilee), roughly 120 km north to the Dead Sea along the Jordan Valley and 


212m BSL (Alpert et al., 1982). In this study, model simulation of the air flow in Lake 
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Kinneret identified the strong entrance of the Mediterranean Sea breeze to the western 
shore of the lake, while on the eastern shore the winds were much weaker due to the 
unique topography of the area. This effect is caused by the interaction of the synoptic 
system with the meso-scale sea-breezes. The significant weakening on the eastern side 


does not evoke the same kind of emissions as showed by Kishcha et al. (2013). 


The use of normalized AOD difference (NAODD index) made it possible to 
capture the seasonal variability in the aerosol levels between the eastern and western 
coasts of the Dead Sea. Mishra et al. (2014) studied aerosol properties over the eastern 
Mediterranean basin using a verity of techniques. They found that the summer 2010 AOD 
means are 0.22 + 0.02 for polluted dust, 0.11 + 0.04 for polluted continental, 0.10 + 0.04 
for dust, and 0.06 + 0.01 for marine aerosol. In comparison, our results show that for the 
summer of 2013, the AOD mean was 0.174 0.01 suggesting of polluted dust. In the 
future, adding ground monitoring stations would help to further validate our analyses. 
MODIS 3 km product should also be included for analyses and compared by means of spatial 
variability in AOD levels. Furthermore, additional research is required to improve the 
accuracy of the AOD estimates above desert regions such as the Middle East (Sayers et 
al., 2014 and De Meij et al., 2011) and in other regions such as South and Southeast Asia 


(Kanniah et al., 2016). 


Since our analysis pertained to one year, more data is required in order to fully 
understand the aerosol regime in the region, the prevailing aerosol patterns, and the 
atmospheric conditions that contribute to high aerosol concentrations. In addition, further 
research would be required in order to better understand the impact of the local industrial 


sources and their dynamics. In our study, we found that there is a “hot spot” area of 
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higher AOD values surrounding the industrial region of Rotem factory (14 pixels, marked 
on Figure 2), reaching a difference of 13% during the winter. Several studies examined 
the impact of local emissions in the New England Region of US (Tang et al., 2017; 
Schutgens et al.,2012). Unfortunately, no studies (to our knowledge) examined the 
impact of local emissions or other anthropogenic activities, e.g. the role of land 
cultivation, population density, human activities and vehicle traffic, in the DS region. In 
this regard, MAIAC represents a strong data-base source for such endeavors. First, its 
high spatial coverage enables data analysis in areas with limited monitoring stations at 
ground level, e.g. the eastern shore of the Dead Sea. Second, its better resolution (1x1 km 
pixels) can identify local sources of pollution. Yet, this may mostly apply for cases with 
high and strong pollution emissions at the ground, which otherwise would have been 
significantly downgraded or gone unnoticed in lower resolution retrievals (i.e. MODIS 
standard10x10 km pixels). Chudnovsky et al. (2013) have shown how with the help of 
MAIAC, spatial variability in pollution levels over the bright urban area of Boston was 


recognized. 


To our knowledge, this is the first study that was able to show AOD patterns over 
the Dead Sea with such clarity, thanks to the MAIACs high 1km spatial resolution. Also, 
at present, 1km model results to the region are not yet available. Kishcha et al., (2016) 
have conducted model simulations in the region for the March 2013 dust event, using the 
online-coupled weather forecast model COSMO-ART, with a grid of 0.025° (~3 km), 
identifying that the maximum near surface dust concentration was located on the western 
part of the valley. Yet this is a description of a singular event (‘case study’), similar to the 


work of Gasch et al. (2017) (a grid of 2.5km) and Alpert et al. (2016) (a ceilometer study) 
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focusing on the September 2015 dust event. While our work has shown such a trend 
following a year of data analysis, which in turn allows us to classify climatological 


patterns. 


Using both passive and active means of remote sensing at a high spatial and 
spectral resolution with angular observations above desert regions is the next step in 
improving current knowledge about dust and anthropogenic pollution. In this regard, the 
integration of a variety of satellite data on a pixel basis - combining data from passive 
systems such as OMI, MODIS, and VIIRS (which provide the spatial pattern of 
pollution) with data from active systems like CALIPSO (which provide information 
about the vertical distribution and classification of aerosols) - is required, given the close- 
knit relation between aerosol levels and atmospheric circulation, pollution sources, and 
radiation (e.g. Gkikas et al., 2016; Pappalardo et al., 2014; Thies and Bendix 2011). 
Establishment of additional ground measurement sources such as ceilometers and Lidars 


will help to fully analyze the 3D aerosol regime in the DS area. 


5. Concluding remarks 


Until recently, the aerosol patterns above bright desert regions at high spatial 
resolution were unavailable. In this paper we used | km AOD retrieval from MODIS data 
based on the MAJAC algorithm and generated for the first time the prevailing aerosol 
spatial patterns above the Dead Sea coasts on a seasonal basis. These patterns are 
governed by the interaction of the synoptic system with the meso-scale sea-breezes and 
the region complex topography. Furthermore, we have shown that the high spatial 
resolution is essential to identify industrial zones/ emissions as apparent from our 


preliminary analysis over Rotem factory. 
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Despite promising results, more data need to be processed and analyzed to 
understand the full potential and limitations of the high resolution MAIAC AOD product 
over bright desert regions. This will allow us to improve the accuracy in AOD retrieval 
and ground pollution level estimations and will greatly assist in future policy to preserve 


the unique Dead Sea area. 
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Tables 


Table 1. Seasonal comparison between AERONET data and MODIS data (DT, DB, and 
MAIAC) for 2013 for the three AERONET satiations in Israel: Eilat, Sde Boker, and Nes 
Ziona. 


Nes Ziona Comparison - AERONET vs MODIS (DT, 


DB, MAIAC) 
Data source Statistics | Fall | Winter | Spring | Summer | Yearly 
DT N 28 22 34 48 132 
R 0.81 | 0.81 0.88 0.56 0.80 
Intercept | 0.12 | 0.05 0.08 0.15 0.10 
Slope 0.56 | 1.10 1.08 0.74 0.91 
DB N 36 34 47 68 186 
R 0.53 | 0.68 0.86 0.39 0.75 
Intercept | 0.07 | 0.07 0.06 0.11 0.02 
Slope 0.39 | 0.89 1.06 0.44 0.86 
MAIAC N 31 24 38 75 168 
R 0.80 | 0.77 0.6 0.53 0.64 
Intercept | 0.02 | 0.07 0.08 0.04 0.06 
Slope 0.88 | 0.66 0.52 0.73 0.68 


Sde Boker Comparison - AERONET vs MODIS (DB, 


MAIAC) 
Data source Statistics | Fall | Winter | Spring | Summer | Yearly 
DB N 67 41 60 80 248 
R 0.26 | 0.27 0.62 0.27 0.51 
Intercept | 0.23 | 0.21 0.30 0.38 0.28 
Slope 0.36 | 0.37 0.46 0.32 0.26 
MAIAC N 56 35 51 79 221 
R 0.64 | 0.60 0.57 0.79 0.60 
Intercept | 0.03 | 0.03 0.11 0.03 0.09 
Slope 0.65 | 0.58 0.27 0.61 0.33 


Eilat Comparison - AERONET vs MODIS (DB, MAIAC) 


Data source Statistics | Fall | Winter | Spring | Summer | Yearly 
DB N 26 14 53 52 145 
R 0.57 | 0.97 0.59 0.62 0.48 
Intercept | 0.07 | 0.03 0.11 0.03 0.12 
Slope 0.25 | 0.53 0.24 1.18 0.29 
MAIAC N 23 11 50 58 142 
R 0.68 | -0.15 0.65 0.66 0.58 
Intercept | 0.11 | 0.17 0.19 0.11 0.17 
Slope 0.6 | -0.11 0.20 0.76 0.25 


Table 2. Seasonal comparison between DT/DB and MATAC AOD for the same days 


and locations. 


Data source Season Variable DT DB MAIAC 
Nes Ziona R 0.57 0.47 0.59 
AERONET neN=2) p-value 0.0016 0.047 0.0017 
; R 0.478 0.439 0.62 
Spring (N=26) | vatue | 0.0009 | 0.0028 0.0009 
R 0.48 0.32 0.26 
pee). | aaavatie. |, 20/006 0.0034 0.25 
R 0.94 0.87 0.88 
Winter (N=22) |  aiue | p<0.0001 | p<0.0001 | 0.0001 
Sde Boker Fall (N=61) R ---- 0.24 0.63 
AERONET p-value ---- 0.063 p<0.0001 
Spring (N=53) R ---- 0.68 0.72 
p-value ---- p<0.0001 | p<0.0001 
Summer (N=84) R no-- 0.26 0.70 
p-value ---- p<0.024 p<0.0001 
Winter (N=36) R ---- 0.27 0.56 
p-value ---- 0.10 p<0.0001 
Eilat Fall (N=34) R ---- 0.29 0.75 
AERONET p-value ---- 0.09 p<0.0001 
Spring (N=48) R ---- 0.47 0.67 
p-value ---- 0.0007 p<0.0001 
Summer (N=51) R ao-+ 0.46 0.63 
p-value ---- 0.0007 p<0.0001 
Winter (N=10) R ---- 0.11 -0.0059 
p-value ---- 0.78 0.98 


Table 3. Comparison of pollution levels between the eastern and western coasts of the Dead 
Sea for 2013, by season. Pollution focal, indicating the coast with the highest pollution on a 
given day, was determined by using the Normalized AOD difference index (NAODD). 
NAODD positive stands for higher aerosol levels above the east, while NAODD negative 
indicates higher aerosol levels over the western coast. 


Pollution focal: 
Season NAODD Frequency Percent (%) 
NAOOD>0O (East) 29 37.66 
Fall (77 days) 
NAOOD<0O (West) 48 62.34 
NAOOD>0 (East) 29 43.28 
Spring (67 days) 
NAOOD<O (West) 38 56.72 
NAOOD>0O (East) 7 7.69 
Summer (91 days) 
NAOOD<0O (West) 84 92.31 
NAOOD>0 (East) 24 53.33 
Winter (45 days) 
NAOOD<0O (West) 21 46.67 
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Figure 1: AOD for selected days in August above the Dead Sea region. Panel A shows RGB 
data; Panel B shows MAIAC AOD data at a 1km resolution; and Panel C shows MODIS 


Collection 6 combined (Dark Target + Deep Blue) AOD product. Note missing AOD values 


when compared to MAIAC. 
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Figure 2: Study Area. Left: a map of the area of Israel, with the Dead Sea region highlighted 
by the dotted line (major cities, AERONET stations, and Masada historical site marked); 
Right: an enlargement of the study area (rural settlements, cities, and industrial zones 


marked). The historical site of Masada is also marked. 
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Correlations were computed between 
MAIAC, MODIS DT and MODIS DB 
retrievals to AERONET data to compare 
their strength of retrieval over the region 
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Case Study: Analysis Comparison between Calculation of NAODD index for Production of AOD 
of the industrial zone east and west coasts the year 2013 to identify east/west maps for typical 
Rotem Factory of the Dead Sea differences between the 2 shores aerosol patterns 


Figure 3: Flow chart of the data analysis process. 
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Figure 4. Monthly trend of AOD for Sde Boker: Comparison between AERONET and 


MAIAC AOD retrievals. Error bars represent standard deviation. 
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Figure 5. Seasonal average of MAIAC AOD for 2013 (where N represents the number of 


available days for analysis in each season). 
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Figure 6: AOD for selected days in above the Dead Sea region, each one representative of a 
different typical AOD pattern: (A) A clear day with hardly any aerosols over the region, on 
July 2, 2013 (the dotted line represents the division between the eastern and western parts of 
the Dead Sea region); (B) High AOD values concentrated along the west coast only (on July 
4, 2013); (C) High AOD on the eastern side of the Dead Sea (March 7, 2016); and (D) 
elevated levels of aerosols above the entire region, with significantly higher aerosols 
concentration surrounding both eastern and western coasts of the Dead Sea (on the September 


17, 2013). 


Highlights 


e Relatively good performance of MODIS MAIJAC AOD retrievals over arid 
zones. 
Different AOD patterns emerged over the Dead Sea coasts. 
Western and eastern coasts showed seasonal dependence in aerosol 
concentrations. 

e ‘Rotem’ industrial zone was clearly identified. 


